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An accurate quantification and characterization of the available wind resources is necessary to 
optimally design a wind farm. To effectively evaluate the wind energy, studying the wind’s statistical 
characteristics is required. The probability distribution of wind speed is a very important piece of 
information needed in the assessment of wind energy potential since wind power is proportional to the 
cube of wind speed. Therefore, choosing a probability function having high goodness of fit with the 
observation data plays a quite significant role in wind energy assessment. In this study, three 
probability density functions, i.e., two-parameter Weibull, Logistic and Lognormal are employed to 
wind speed distribution modeling using data measured at a typical site in Inner Mongolia, China, over 
the latest three year period from 2009 to 2011. The performance of these three functions is compared 
so as to select the best one. As one of the most favorable distributions, Weibull function is a most 
applicable approach in describing wind speed’s distribution in many cases. However, though 
performances of three particle swarm optimization algorithms and 18 differential evolution approaches 
of estimating the shape parameter estimation in the Weibull function are compared, then the one 
which performs best is selected to determine the optimal shape parameter to obtain the most accurate 
shape parameter estimation results. The performance of the Weibull function is worse than the Logistic 
under the measured wind speed data and the chosen error evaluation criteria. Besides, as compared to 
the Lognormal function, the Logistic function provides a more adequate result in wind speed 
distribution modeling. Therefore, in this work, the Logistic function is applied to the consequent wind 
energy assessment through the availability factor, capacity factor, and turbine efficiency of a wind 
turbine. Assessment results have shown that it is suitable to build a wind farm in this area. 

© 2013 Elsevier Ltd. All rights reserved. 
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1. Introduction 


With the fast development of China’s economy over the last 
few decades, China is becoming a big power in energy consump- 
tion. The consumption of coal, petroleum and nature gas con- 
stituted 74%, 18.8%, and 4.1% of the total annual primary energy 
consumption in 2009 [1], respectively, and the corresponding 
proportions have changed to 68%, 19% and 4.4% in 2010 [2]. 
Though the total proportion of these three fossil fuels decreased 
from 96.9% in 2009 to 91.4% in 2010, the total annual primary 
energy consumption has increased from 2920.28 million tons of 
standard coal as calorific value calculation in 2009 to 3249.39 
million tons in 2010, which means the consumption of these 
three fossil fuels in 2010 is still larger than that in 2009. As 
known, fossil fuels are non-renewable resources. Meanwhile, a 
large amount of greenhouse gases is exhausted in the burning of 
the fossil fuels, which causes an increasing crisis in global climate 
warming. Thus, in order to deal with the energy crisis and global 
climate warming, clean and renewable energy is urgently needed. 
The “Medium and long term plan of the renewable energy of 
China” proposed that it is our goal to increase the renewable 
energy consumption of China up to 15% in 2020 [3]. As one of the 
cleanest and renewable energy resources, wind energy has been 
drawing the attention of people, and a large number of wind 
farms have been built for wind power utilization. Meanwhile, 
new areas to design wind farms should still be explored. To 
optimally design a wind farm, an accurate wind energy resource 
assessment is necessary. It is not only an essential part of the 
development of wind power utilization, but also provides the 
investors with the necessary confidence in financial feasibility 
and mitigating risks [4]. 

To effectively evaluate the wind power available for a parti- 
cular site, studying the wind’s statistical characteristics is needed. 
For this reason, a large number of studies have been published 
concerning the use of a variety of probability density functions 
(pdfs) to describe wind speed frequency distributions. For exam- 
ple, five distributions including the Weibull, Rayleigh, gamma, 
lognormal, and inverse Gaussian pdfs were used by Zhou et al. [5] 
to investigate the wind speed distribution for wind speed data 
measured from five representative sites in North Dakota. Another 
two new pdfs named Pearson type V and Burr were introduced by 
Valerio et al. [6] to describe wind speed frequency distributions in 
the urban area of Palermo, in the south of Italy. Different from 
them, the performance of another pdf: the binormal distribution 
function was tested by Peter and Imre [7] through using ERA-40 
surface wind speed data over Europe. Considering that the single 
Weibull function may not best reflect the distribution of the wind 
speed, a bimodal Weibull & Weibull pdf was developed by 


Jaramillo and Borja [8] to analyze the wind speed frequency 
distribution in La Ventosa, Mexico. In addition, two researchers 
[9] developed a theoretical approach to the analytical determina- 
tion of wind speed distributions of four stations in Elazig through 
application of the maximum entropy principle in 2007, and in 
2009, the same two authors [10] represented another theoretical 
approach of wind speed frequency distributions through applica- 
tions of a Singly Truncated from below Normal Weibull mixture 
distribution and a two component mixture Weibull distribution. 
Furthermore, the minimum cross entropy (MinxEnt) probability 
density function (pdf) derived from the MinxEnt principle, which 
covers the maximum entropy principle, was used to determine 
the diurnal, monthly, seasonal and annual wind speed distribu- 
tions by Kantar and Usta [11]. 

Besides, for a specific pdf, methods adopted to determine the 
parameters in this pdf are various. As one of the most popularly 
used models due to its ability to fit most accurately the variety of 
wind speed data measured at different geographical locations in 
the world [9], the Weibull function is frequently used to the wind 
speed distribution modeling. Methods applied to Weibull para- 
meters estimation commonly including: the graphical method 
[12], the maximum likelihood method [13], the moment estima- 
tion and quartiles [14], as well as the modified maximum like- 
lihood and energy pattern factor method [15], etc. 

In this paper, three probability density functions including the 
Weibull, Logistic and Lognormal functions are employed to wind 
speed distribution modeling. To obtain the most accurate shape 
parameter value in Weibull function, three particle swarm opti- 
mization algorithms and 18 differential evolution approaches are 
adopted, and the one with the best performance is selected to 
determine the final optimal shape parameter. Then the perfor- 
mance of these three probability density functions is evaluated 
according to four error evaluation criteria by giving them differ- 
ent weight, the optimal function in modeling the wind speed 
distribution is chosen to assess the wind energy of this area. 

The rest of this paper is organized as follows: Section 2 
presents the wind energy potential assessment review. In 
Section 3, the extrapolated data in the hub height are presented 
as well as a brief analysis to the data. Theory of the two- 
parameter Weibull function is introduced in Section 4, and the 
objective function which will be used in the subsequent para- 
meter optimization process is also exhibited. Then the algorithms 
of three particle swarm optimization and 18 differential evolution 
approaches are presented in Section 5. Section 6 gives the final 
optimal results in estimating the shape parameter in the Weibull 
function. The performance of two other probability density 
functions are compared with the Weibull function under four 
error evaluation criteria in Section 7, and the best function 
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selected in Section 6 is applied to the wind energy assessment in 
Section 8. Finally, Section 9 concludes the findings of this paper. 


2. Wind energy potential assessment review 
2.1. Assessment related to the Weibull distribution 


Due to the superiority performance of the Weibull distribution 
in describing the wind speed distribution, many wind energy 
potential assessment works have been carried out with the help 
of the Weibull distribution. 

Among these works, periods of the wind speed data covering 
are various. Weisser [16] utilized the Weibull density function to 
estimate the wind energy potential in Grenada, West Indies by 
the wind speed data from 1996 to 1997. A data source consisting 
of five years (1998-2002) of hourly mean wind data was adopted 
to determine the wind energy potential in Maden-Elazig, Turkey 
[17] with the assistance of the Weibull distribution. The same 
distribution was also used by Keyhani et al. [18] to find out the 
wind energy potential of the capital of Iran, Tehran according to 
the statistical data of eleven years’ wind speed measurements. 
Ohunakin and Akinnawonu [19] examined the wind energy 
potential in Jos, Nigeria using the wind speed data during 1971- 
2007 subjected to the Weibull analysis as well. Moreover, Chang 
et al. [20] analyzed the wind characteristics in Taiwan based on 
an even long period from 1961 to 1999 so as to present an 
assessment of the wind energy potential in this area. 

In addition, the characteristics of the wind speed data at 
different heights can also be described by the Weibull distribution 
so as to provide a reliable wind energy potential assessment. The 
possible wind energy potential of Kütahya, Turkey at the single 
30m height was investigated by Köse [21] according to the 
Weibull distribution. Ullah et al. [22] employed the Weibull 
distribution to provide an evaluation of wind energy potential 
in Kati Bandar, Pakistan at three different heights, i.e., 10 m, 30 m 
and 50m. The Weibull distribution was also employed by 
Mirhosseini et al. [23] to assess the wind energy potential 
locations in province of Semnan in Iran with the wind speed data 
at 10m, 30m and 40m height. And based on the wind data 
collected at the heights of 10, 40, 50 and 60m, Đurišić and 
Mikulovi¢ [24] carried out an assessment on wind energy in the 
South Banat region, Serbia. 

As a special form of the Weibull function, the Rayleigh 
distribution function has also been frequently used to assess the 
wind energy potential. A Rayleigh distribution was adopted for 
defining the distribution of wind velocity in terms of its prob- 
ability density and cumulative distribution functions by Mathew 
et al. [25] so as to present a wind energy potential assessment of 
eight sites in Kerala, India. Rayleigh parameter at a height of 30 m 
has been estimated and used by Amar et al. [26] to assess the 
wind energy potential of Sidi Daoud, Tunisia. 


2.2. Assessment related to other approaches 


Apart from the probability density functions, some other 
approaches have been applied to the wind energy potential 
assessment as well. Onat and Ersoz [27] analyzed the wind energy 
potential of three regions in Turkey with the help of the Wind 
Atlas Analysis and Application Program software, the same soft- 
ware was also used by Durak and Sen [28] to assess the wind 
power potential in Turkey and Akhisar. Hossain et al. [29] 
assessed the wind energy potential for India by an innovative 
approach using the geographical information system (GIS) plat- 
form. A wind resource inventory created by the GIS data supplied 
with the wind maps which was developed with a meso-scale 


wind energy model coupled with wind data was presented by 
Acker et al. [30] to assess the wind energy potential in the state of 
Arizona as well. And Thiaw et al. [31] used the multi-layer 
perceptron type artificial neural networks to determine the wind 
energy potential on the Dakar site. Besides, the fifth-generation 
meso-scale Model was used in conjunction with an h-adaptive 
finite element model by Pepper and Wang [32] to assess the wind 
energy potential of Nevada, and an assessment of long-term wind 
energy potential at the demonstration offshore wind farm in 
Korea was estimated by Oh et al. [33] by using the Measure- 
Correlate—Predict technique. 


3. Site representation and brief analysis 
3.1. Site representation 


Wind speed data used in this paper is measured from 2009 to 
2011 at one of the sites in Inner Mongolia, China, with the latitude 
of 42.5°N and longitude of 115°E. The data are sampled at the 
standard height 10 m with a temporal resolution of 6 h, so there 
are four data for one day. The cut-in, rated and cut-off speeds of 
the wind turbine S 43/600 in this area are 3, 15 and 24 m/s, 
respectively, and the rated power is 600 kW. Rotor’s diameter is 
43 m, which has a blade sweep area of 1452.2 m°, and the hub 
height is 50 m above the ground level. 

Wind speed values sampled at the standard height are first 
extrapolated to the ones at the hub height using the following 
empirical power law: 


v h\” 

n= (70) z 
where v and v are the wind speeds at the standard height 10 m 
and at the hub height h m, respectively, and « is the roughness 
factor, in this analysis, a typical value of 0.14 is adopted. By 
substituting h=50 and the measured wind speed data at the 
standard height 10m to Eq. (1), wind speed values at the hub 
height are obtained. Fig. 1 shows the variation trend of the wind 
speed throughout the year. 

As demonstrated in Fig. 1, the highest and lowest wind speeds 
in this area are about 18 and 2 m/s, respectively, which appear in 
February and May. In addition, except for May, June is also one of 
the low wind speed period. 


T T 
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Fig. 1. The variation trend of the wind speed throughout the year. 
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Table 1 
Brief analysis of the wind speed characteristics. 


Characteristic Period 

January February March April May June July August September October November December Yearly 
Mean speed (m/s) 8.9439 10.2558 8.6957 8.2298 7.4566 7.5508 7.8993 9.6786 9.5821 9.5605 9.6972 8.7084 8.8444 
Standard deviation (m/s) 1.9704 2.6717 1.8471 2.0910 1.8862 1.9052 2.0251 2.3573 2.3878 2.5014 2.6589 2.0899 2.3738 
Skewness —0.512 0.374 0.198 0.670 -—0.288 0.712 0.027 -—0.071 0.178 0.373 —0.049 0.196 0.338 
Kurtosis 0.526 -—0248 -—0.188 0.530 -—0.252 1.144 -0407 -—0.379 0.188 —0.361 —-—0.674 0.003 0.067 


3.2. Brief analysis to the wind speed data 


In this section, the wind speed data at the hub height are 
analyzed from the following three aspects: the central tendency, 
the discrete degree and the distributional pattern. 


3.2.1. Introduction to the statistics 


3.2.1.1. Statistics used to describe the central tendency. The central 
tendency refers to the tendency of a set of data which draws 
close to a central value. Statistics used to describe the central 
tendency reflects the representative value or the center value of a 
data sequence. The mean value of a data sequence is usually 
employed to describe the central tendency, which is defined as 
follows: 
n 

Xi, (2) 


i=1 


sie 


X= 


where X is the mean value, {x;}7_ , is the data sequence and n is 
the number of the data contained in this data sequence. 


3.2.1.2. Statistics used to describe the discrete degree. Using only 
the mean of a data sequence cannot depict the property of the 
data completely. Apart from the central value, the ‘distance’ of all 
data relative to the central value should also be investigated. If 
the data closely gather around the central value, the discrete 
degree of the data is small, which means the central value is the 
representative of the data sequence; Otherwise, if the data are 
loosely distributed around the central value, the discrete degree 
of the data is large and the central value cannot represent all of 
the data well in this case. Therefore, the central value and the 
discrete degree of the data should be given together to give a 
comprehensive description to the data sequence. The discrete 
degree of the data sequence relative to the mean value can be 
described by the Standard deviation (Std) value expressed as 
follows: 


Std = (3) 


where x is defined by Eq. (2). 


3.2.1.3. Statistics used to describe the distributional pattern. To 
grasp the characteristics of the data better, the distributional 
pattern of the data should also be surveyed. The distributional 
pattern can be researched by two statistics: Skewness and Kurtosis, 
the expression of these two statistics are shown in the following 
equations: 


1< E E 
Skewness = repr (x;-xy /Std?, (4) 


i=1 


: 1 = fej a4 
Kurtosis = a X @%-%)*/Std*-3, (5) 


i=1 


where X and Std are defined by Eqs. (2) and (3), respectively. 
Skewness is a description of the symmetrical characteristic of the 
data, where Skewness=0 means the distributions of the data is 
symmetric, or else, it is not symmetric, and the larger the absolute 
of the Skewness value, the more skewness of the data. Kurtosis is a 
statistic employed to describe the steep degree of the data, where 
Kurtosis=O representing the steep degree of the data is the same 
as the standard normal distribution, Kurtosis >0 means the 
distribution of the data is steeper than the standard normal 
distribution while Kurtosis<0O reflects the distribution of the 
data is less steep than the standard normal distribution. 


3.2.2. Analysis results to the wind speed data 

Applying the above statistics introduced in Section 3.2.1 to the 
wind speed data, results listed in Table 1 are obtained. 

As demonstrated in Table 1, the mean speed in February is the 
maximum, while the one in May is the minimum. Wind speed 
data in November present the highest degree far away from the 
mean value for the largest standard deviation value it reached, 
and the discrete degree of the wind speed in March is the lowest. 
Besides, the skewness of the wind speed data in June is the most 
severe, and wind speeds in January, April, June, September, 
December and the whole year are steeper than the standard 
normal distribution (for the Kurtosis values of these months are 
larger than 0), while the wind speed data in other months are less 
steep than the standard normal distribution. 

Surveying the wind speed characteristics and distribution 
according to these statistics is not enough. More detailed infor- 
mation should be acquired by analyzing the wind speed distribu- 
tion and fitting it through some distribution functions, such as the 
two-parameter Weibull model, the Logistic and the Lognormal 
distributions, etc. 


4. Two-parameter Weibull model 
4.1. Weibull function 
Weibull probability density function (pdf) is expressed as 


w- (GO 0] © 


where v, k, and c are the wind speed, the dimensionless shape 
parameter, and the scale parameter with the same unit as the 
wind speed, respectively. Based on the Weibull pdf, Weibull 
cumulative distribution function (cdf) can be calculated by the 
integration: 


ror [soe OO O] 
=- | dexp |- È) i = —exp |- (-) i 
=1-—exp |- ©) j ; (7) 


J. Wu et al. / Renewable and Sustainable Energy Reviews 21 (2013) 215-228 219 
Table 2 
Simulation results by the three PSO algorithms. 
Algorithm Item Minimal iteration when the termination criterion is the fitness level is 1 x 107" 
January February March April May June July August September. October November December Yearly 
BPSO Iteration 29 28 50 43 31 54 88 38 50 42 43 42 50 
k 5.2172 4.3412 5.4290 4.4622 4.4836 4.4965 4.4184 4.6744 4.5584 4.3205 4.1024 4.7513 4.2007 
Fitness (x 107°) 6.0542 6.3340 7.1916 1.7391 3.6284 1.0514 1.8565 3.0458 1.8726 5.0785 4.5689 9.2880 8.1431 
WPSO Iteration 21 24 26 17 19 18 22 17 20 16 14 24 15 
k 5.2175 4.3408 5.4290 4.4622 4.4835 4.4961 4.4180 4.6749 4.5583 4.32 4.1027 4.7505 4.2006 
Fitness (x 107°) 1.7934 3.7623 6.5053 3.3368 5.4070 9.6330 7.2239 9.3189 0.1790 7.7823 3.7711 8.7104 4.5387 
QPSO a _a a a a a a a a a a a a a 
Fitness value when the termination criterion is the maximum iteration number equals to 50 
BPSO k 5.2174 4.3408 5.4296 4.4620 4.4838 4.4964 4.4183 4.6742 4.5585 4.3202 4.1026 4.7505 4.2005 
Fitness (x 107°) 2.8315 2.9134 1.8608 2.5735 0.5046 1.9653 0.2005 7.3031 4.4822 3.3872 0.3628 7.6569 2.0268 
WPSO k 5.2176 4.3409 5.4295 4.4621 4.4838 4.4964 4.4183 4.6745 4.5583 4.3203 4.1025 4.7509 4.2004 
Fitness (x 1078) 0.1586 0.7320 1.5120 1.0433 0.1898 0.0651 0.0439 0.0437 0.4090 0.0703 0.2473 0.0121 0.0403 
QPSO k 5.2176 4.3409 5.4295 4.4621 4.4838 4.4964 4.4183 4.6745 4.5583 4.3203 4.1025 4.7509 4.2004 
Fitness (x 1071°) 0.3651 0.0872 0.0926 0.0401 0.2381 0.0183 0.1014 0.8759 0.0230 0.0348 0.0541 0.0002 1.7107 
è Incalculable value. 
4.2. Objective function 5. Parameter optimization approaches 


Based on the Weibull pdf, the mean 7 and standard deviation 
(o) can be obtained: 


T= m uf (v)dv 


L O0 e-@ | 


= [exe |-) Je) 
oT ORA 


ep (1 + z) (8) 
o? = Ev—v)" = Ev? —(Ev)? 

-[{- 2) @° |-@ o-er) 

=f ef w O pOr) 

=er(1+3) -er (1+7): (9) 


To search the Weibull parameters by using the intelligent 
parameter optimization approaches, the residual value e defined 
as below is used as the objective function in this paper just as Liu 
et al. did in [34]: 


o2 T(1+2/k)-r°(1+1/k) 
P T? (1+1/k) 


; (10) 


If the objective function converges to a predefined tolerance 
value then an acceptable shape parameter k is obtained; further- 
more, the scale parameter c can be determined using Eq. (8). In 
the process of the parameter searching, when a given maximum 
iteration number has been produced or a satisfactory fitness level 
is reached, the final value of the parameter is determined and the 
algorithm is terminated. 


To search the optimal shape parameter in the Weibull func- 
tion, the following intelligent parameter optimization approaches 
are employed. 

5.1. Three particle swarm optimization (PSO) algorithms 
5.1.1. Basic PSO (BPSO) 

PSO is a population based optimization tool, where the 
system is initialized with a population of random particles and 
the algorithm searches for optima by updating generations. 
Suppose that the search space is n-dimensional, the ith particle 
is represented by an n-dimensional vector X;=(Xj1,Xi2,...,Xin) and 
the velocity of this particle is represented as V;=(vj1,Vj2,...,Vin). 
According to the objective function, the fitness of each particle 
can be evaluated. Record the fitness of the personal best position 
P;=(Di1,Di2,...,Din), Which is the best previously visited position of 
the particle i, as fpbes- Meanwhile, the position of the best 
individual of the whole swarm is marked as the global best 
position Py=(Pg1,Pg2,...,Psn), and the objective function 
corresponding to Pg is fobesr. At each step, the velocity and 
position of the ith particle are updated according to the 
following two equations until a user-defined stopping limit is 
reached: 


V; = ORV; + C1 71 *(P; —X;) + C2*r2*(G—X;), (11) 


X;=X;4+ Vi, (12) 


where œ is called the inertia weight, by which the impact of 
previous velocity of particle on its current one can be controlled, 
rı and rz are independently uniformly distributed random variables 
with range (0, 1) and cı and c are positive constant parameters 
called acceleration coefficients which control the maximum step 
size. In the BPSO algorithm of this paper, œ is a random value in the 
range (0, 1). 


5.1.2. PSO with self-adaptive stochastic inertia weight (WPSO) 

As exhibited in Eq. (11), the velocity of particle is influenced by 
the inertia weight œ. Different from the BPSO, in the WPSO 
algorithm, the inertia weight œ is defined as values in different 
ranges according to the change of fopest: 


(a) If fobest is unchanged, œ is defined as a random value in the 
range (0.5, 1) to increase the search intensity. 
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(b) Otherwise, if fgbest is changed, œ is defined as the same as in 
the BPSO algorithm, i.e., it is defined as a random value in the 
range (0, 1) [35]. 


5.1.3. Quantum-behaved particle swarm optimization (QPSO) 

QPSO algorithm is a new global optimization algorithm based 
on the PSO algorithm. Different from the PSO algorithm, only the 
location information for the particles is given in QPSO algorithm, 
the location is updated by the following three equations: 


Pia = P*pbestia +(1—p)*gbesta, (13) 
1 

mbest = M> pbest;, (14) 

Xid = Dig + 08| mbesty—Xia|*In(1/u), (15) 


where ġ and u are random values in the range (0,1), thus pig is a 
random position between pbest and gbest, M is the size of the 
particles, and «œ is called the expansion-contraction factor. As 
shown in Eq. (14), mbest is the center of the current optimal 
location in all individual particles. Parameter « in Eq. (15) is used 
to control the convergence speed of the algorithm. It is the sole 
parameter needed to determine and control, usually it is deter- 
mined as the linear function of the iterations, that is 


& = (max —%min)**(iteration—iter) /iteration + &pin (16) 


where iter is the current iterations, iteration is the total iterations, 
Xmax and Xmin are two positive constants, which are usually set as 
1.0 and 0.5 [36]. 


5.2. Differential evolution 


In the algorithm of the differential evolution (DE), a set of D 
optimization parameters is called an individual, which is repre- 
sented by a D-dimensional parameter vector, and the population 
consists of NP parameter vectors X¢ (i=1,2,...,NP for each gen- 
eration G). The basic process of DE includes three steps: mutation, 
crossover and selection, which are described as follows: 


Step 1, Mutation: For each target vector X¢ (i=1,2,...,NP), the 
mutation operation is conducted according to Eq. (17) so as to 
generate a mutant vector V¢*!: 


vot! — x6 4 Fx(xS xG), r1 412 413 Ai, (17) 


where r1, r2, r3 are integer indexes randomly chosen from 
{1,2,...,NP}, F is called the mutation factor and is taken from 
the interval [0,1]. Actually, it is a factor which controls the 
amplification of the difference between two individuals. In this 
paper, this mutation method is called Strategy 1. In addition, 
eight other mutation approaches which are represented as 
follows are employed to optimal parameter searching in this 
paper: 

Strategy 2: v+! = xf + Fy (x6, xf) + Fox (xS -x5 ); 

Strategy 3: v+! =x + (x —xS )*((1 -0.9999)*rand +F); 

Strategy 4: vS +! = x$ + Fy(x5—x&), Fi. =(1 —F)*rand+F, and 

F; are equal for all of the parameters needed to estimate; 

Strategy 5: v+! =x +Fyx(xG—xS), F)=(1—F)*rand+F; 

Strategy 6: v+! = xf + F*(x6, -x6 ); 

Strategy 7: vot! = xG + Fae(xS -x6 +x9,—xS): 

Strategy 8: vO+! = xg + F(xS -x6 +x9 -x5 ); 


Strategy 9: 
an xG +F*(x6 -x6 ) if rand < 0.5, 
i XE +0.5x(F + 1)(xG +xG—2%xS) if rand > 0.5, 


where r4 and r5 are integer indexes randomly chosen 
from 1,2,...,.NP as well, x¥ and xf... are the ith and the best 
individuals in generation G, respectively. 
Step 2, Crossover: In order to increase the diversity of the 
perturbed parameter vectors, crossover is introduced. The 
target vector is mixed with the mutated vector, using the 
following scheme, to yield the trial vector 


U+! = (uf uSt1,...ugt!), that is 


— vat? if rand(j) < CR or j=rnb(i) © r : 
ee xi otherwise p= Dene A9 
where rand(j) is the jth evaluation of a uniform random 
number generator between [0,1]. CR is the crossover constant 
between [0,1] which has to be determined by the user. In Eq. 
(18), rnb(i) is a randomly chosen index from 1,2,...,D which 
ensures that Uf*' gets at least one parameter from V/’*'. 
Otherwise, no new parent vector would be produced and the 
population would not alter. This is called the binomial cross- 
over. Besides this, another crossover approach named the 
exponential crossover is also a commonly used crossover 
method, in which the trial vector US*! is described as 


A vit? ifje{k,<k+1 >n.. <k+L-1 >n} | 
Uji = e 
a9) 


]J= 1,2,. ” aD, 

where ke{1,2,...,n} is a random index, L is a random value in 
{1,2,....n} and <j>nisjifj<nandj—nifj>n. 
Step 3, Selection: Assuming that it is the goal of the objective 
function to obtain the minimal value of one variable, the 
selection operation is carried out according to the following 
rule: 


xer = uy if f(Ur*') <f(xF). (20) 


x¢ otherwise, 


otherwise, 


That is, if vector Uf*! yields a smaller fitness value than Xf, 
then Xf*! is set to USt'; otherwise, the old value Xf is 
retained. The purpose of this step is to guarantee all the 
individuals of the next generation are as good as or better 
than their counterparts in the current generation. 


6. Shape parameter searching results 
6.1. Simulation results by the above three PSO algorithms 


Based on the three PSO algorithms introduced in Section 5.1, 
the shape parameter (k) is estimated under two termination 
criterions: one is the fitness level equals 1 x 107? and the other 
is the maximum iteration number equals 50. Simulation results 
by the three PSO algorithms are listed in Table 2. 

Note that in the QPSO algorithm, the total iterations, i.e., the 
maximum iteration number, must be given before searching the 
optimal parameters, however, if the termination criterion is the 
fitness level other than the maximum iteration number, this 
algorithm cannot be adopted any more. Thus, in Table 2 corre- 
sponding values are incalculable for QPSO when the fitness level 
(1x 107%) is set as the iteration stopping limit. As shown in 
Table 2, though the optimal shape parameters searched by the 
three PSO algorithm are nearly the same, the WPSO and QPSO 
algorithms both obtain better performance in searching the 
optimal shape parameter (k) than the BPSO: fewer iterations or 
more accurate fitness value. Besides, when the termination 
criterion is the maximum iteration number equals 50, the QPSO 
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obtains the highest fitness accuracy for most of the months, 2 
however, the best performance is not suitable for all of the EJS RARSARSAS 
months, thus the best algorithm should be chosen according to z 
the actual requirements. For example, to receive the best para- £ 
meter optimization performance for December, it is better to E| ARAGA RASS & 
choose the WPSO model to reach the minimum fitness value. 5 
Q 
3 
6.2. Shape parameter searching results by the DE approach 2/S8ESSNSR8RA8 
Y 
In this section, the nine mutation strategies mentioned in $ Ringing em ee B 
Section 5.2.1 and the two crossover approaches introduced in 
Section 5.2.2 are combined to generate 18 different DE methods. 2 
As shown in Eq. (18), the crossover constant CR is determined by P cs en Same, Saag Gaseoi6s 
the users, thus in Table 3 the performances of these methods with A | EAEN E SP NE 
two different crossover constants are surveyed: CR=1 and t 
ap NOnraunnytodt 
CR=0.5. Generally, as seen from Table 3, a smaller crossover FIANNA om 
constant needs larger iteration number to reach the fitness level. =| Hee peneses 
This phenomenon is best explained according to the minimum n a 
iteration number reached by the strategy 7-exp in the case of two 5 | SBALARAILS 
different crossover constants, though their comment goal is to > 
reach the fitness level (1 x 107°). Thus, a bigger crossover con- 2/3 SABAR TESS 
stant can be selected to search the optimal parameter more B qgqeayengesnnr 
quickly. However, though the crossover constant 1 is two times e E 
to 0.5, the corresponding minimum iterations under these two 5 DPaennnr oman 
crossover constants do not have this relationship. z 
Similar to the three PSO algorithms, the performance of these 18 A 
DE approaches when the termination criterion is set as the maximum SI Ew O Oy Gy hb 
iteration number equals 50 for each month and the whole year is ala = 
researched. The results are shown in Table 4. Here, only the results in o q I| naorgenoawn a 
the case of CR=1 are shown. As seen, Strategy 2 and Strategy 3 obtain Se EEI A ERE Ae 
more accurate fitness values than other strategies as well. z = 
5 | MRPs aca 
o 
J a 
6.3. Final optimal shape parameter results y 2 
= g onurtrnner aM 
In this paper, the tolerance value of the objective function is g a A R R 
defined to be 1 x 107° so as to obtain more accurate results than E 8 
it did (1 x 107?) in Ref. [34]. In addition, the optimal parameter in £ 5 
each month and the whole year is selected by comparing the 2 Z| 7PPeS™REORR 
three PSO and the 18 DE algorithms. The final optimal parameter Z E 
for each time stage is selected by comparing the minimum 2 8 EEEE ue cen ET 
iteration number. The approach which obtains the smallest = 5 ve eT T 
iteration number is determined to be the best method and then 3 E 
the shape parameter is estimated by this approach. The shape a z 
parameter in each month and the whole year determined by g a) ae scraas 
which algorithm are plotted in Fig. 2, among which, ‘[a]’, ‘[b]’, ‘[c]’, £ 3 
‘[d]’, ‘[e]’ and ‘[f] mean this optimal shape parameter is searched 5 a | 2 eo oe 
by the Strategy 2-Bin in the case of CR=0.5, Strategy 2-Bin in the = > SQenow 
case of CR=1, Strategy 3-Bin in the case of CR=1, Strategy 1-exp 2 F Ji 
in the case of CR=1, Strategy 9-Bin in the case of CR=1 and E 5 | eras 
Strategy 1-Bin in the case of CR=1, respectively. Scale parameters E Fel E haces, ipo 
determined by the corresponding shape parameters and the is = g 
minimum iteration numbers are shown in Fig. 2 as well. a E enooonn N 
As presented in Fig. 2, the shape parameters searched by the z 5 
Strategy 2-Bin in the case of CR=1 and the Strategy 3-Bin in the E a TLAAHRS Sx 
case of CR=1 are the most, and all of the optimal shape oo x 
parameters are searched by the DE approaches. Besides, all of 2 E 
the shape parameters vary near 5, and the scale parameter line is z 2] aeunraens 
located above the shape parameter line. The minimum iteration E _| 2 
number in searching the optimal shape parameters for the yearly Z l s TE EA 
wind speed data is no more than 5, which is also the smallest in 6 |S AEN EN NEA 
all of the shape parameters searching process. Fig. 3 shows the k ESSS88E888 
complete procedure for wind energy assessment, from which it = =- N m ¥ N 
can be seen that the optimal shape parameter determination is a Gal 5 B B B 8 B B 
necessary part of the best probability density function selection, 3 E E EOE EE E 
so as to provide an accurate wind energy assessment. eS 


Strategy 6 Bin 23 
Strategy 7 Bin 28 
Strategy 8 Bin 25 
Strategy 9 Bin 18 
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Table 4 


Simulation results by DE when the termination criterion is the maximum iteration number equals 50. 


Strategy Fitness value 
January February March April May June July August September October November December Yearly 
Strategy 1 Bin (x 107°) 0.1573 0.0028 0.1070 0.0157 1.2493 0.0374 0.2927 0.1141 2.2147 0.0138 0.1207 0.0503 0.3497 
Exp(x107°) 0.0288 0.7259 0.0281 0.0176 0.5038 0.1216 0.3234 0.7072 0.5759 0.0115 0.1450 0.2090 0.0595 
Strategy 2 Bin(x107'7) 0.8157 0.3215 0.0147 0.5550 0.0582 0.3916 0.0456 0.1642 0.1093 0.1437 0.5338 0.0271 0.2491 
Exp (x107!?) 0.0692 0.3393 0.9403 0.5266 0.1646 7.1883 0.3278 0.0291 0.0159 7.4262 0.4241 1.8338 0.0245 
Strategy 3 Bin(x10~!7) 0.0818 0.1382 0.1490 0.0607 0.0130 0.0096 0.0542 0.1126 0.1200 0.7047 0.1018 0.0163 5.9767 
Exp(x107!7) 1.0958 0.1602 1.6932 0.6109 0.0257 1.5624 0.0738 0.6123 0.1136 0.0173 0.0511 0.0482 0.2150 
Strategy 4 Bin (x 107°) 0.2092 0.0524 0.2356 0.1324 0.4212 0.1525 0.0096 0.2009 3.6355 0.1204 0.7823 0.2187 0.2330 
Exp (x 107° 0.5701 0.1092 0.6607 0.0023 0.2526 0.0122 0.2946 0.0526 0.8931 0.0829 0.4739 0.3023 0.0827 
Strategy 5 Bin ( x 1078) 0.0116 0.0104 0.0697 0.0023 0.0199 2.7678 0.0758 0.0012 0.2862 0.0450 0.0149 0.0423 0.0008 
Exp ( x 1078 0.0073 0.0144 0.0141 0.1230 0.0123 0.0422 0.0015 0.0168 0.0002 0.2261 0.1335 0.0013 0.0004 
Strategy 6 Bin ( x 1077) 0.1539 0.0781 2.2850 5.0708 0.9624 0.0977 0.1522 2.0323 0.8490 1.9373 0.1271 0.0386 1.2156 
Exp (x 1077 1.3920 0.2046 0.0901 0.3051 0.6177 1.8480 0.4300 3.0977 0.5155 0.3120 0.5847 0.0331 0.2196 
Strategy 7 Bin (x 1077) 0.0445 0.0184 0.0898 0.0143 0.0315 0.0147 0.0118 0.0086 0.0169 0.0548 0.0175 0.0048 0.0029 
Exp (x 1077 0.1461 0.0059 1.3403 0.0312 0.1595 0.0340 0.0469 0.0007 0.0195 0.0950 0.0085 0.0438 0.0123 
Strategy 8 Bin (x 107°) 0.1829 0.0450 0.0136 1.1316 0.0631 0.0154 0.0563 3.6629 0.9170 0.0681 0.4209 0.6739 0.0070 
Exp (x 107° 0.2409 0.5098 0.0708 0.0201 0.0026 0.4208 0.5867 0.0429 0.1287 0.1023 0.0531 0.0260 0.0071 
Strategy 9 Bin ( x 107°) 0.7274 0.0992 9.8622 0.0012 0.0265 0.8364 0.1521 0.0202 0.0002 0.8163 0.0671 1.6552 0.0211 
Exp (x 107° 0.5127 0.0245 0.0034 0.1061 0.7093 0.0198 0.4501 0.0067 0.0018 0.4457 0.0968 0.267 0.0125 
15 15 
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Fig. 2. Shape, scale parameter and minimum iteration number during different 
time periods. 


7. Optimal probability density function selection 
7.1. Two other distributions used to compare 


To evaluate the performance of the two-parameter Weibull 
distribution in describing the wind speed distribution, another 
two distributions: Logistic and Lognormal distributions are 
adopted. The pdfs of these two distributions are expressed in 
the following equations, respectively: 


exp[-(v-7)/ò] = (21) 
6{1+exp[—(v—1)/5]} 


fiw) = 


1 [Inwv)—A? 
RO= sg pa? 29? \ (22) 


where 6 = /3a/z is called the Logistic scale parameter, 7 and o 
are estimated in this paper by the moment estimation approach, 
i.e., they are estimated by the mean wind speed and the standard 
deviation of the actual wind speed series. Similarly, 2 and ¢ are 
estimated by the mean and the standard deviation of natural 
logarithm of the observed wind speed. 


Wind energy asse: 
the optimal probabil 


function chosen in Step 3 


Fig. 3. Flowchart of wind energy assessment. 


The cumulative distribution functions of the Logistic and the 
Lognormal distributions can be calculated by the integrations of 
Eqs. (21) and (22), which can be expressed as follows: 


1 1 


P= Ty expl-W—0)/0]_ 1+exp( 7/0) eG 
1 #1 #2 /[* 
F,(v) = 5+ az f exp(—t’)dt, (24) 


where x = In(v)—0/V20. 


7.2. Evaluation criteria 


The performance of the above three distributions is evaluated 
according to four criteria: the first is the mean relative error 
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Fig. 4. Frequency of the yearly wind speed located in each wind speed interval: 
(a) three-dimensional view; (b) x-y view; (c) x-z view; (d) y-z view. 


(MRE), which is given as follows: 


_ 1 z |Xoi—Xci| 
MRE n2 ot (25) 
where {Xoi}; 1 and {x,;}?_ , are the observed and calculated data 
series respectively, and n is the number of the data. Besides this, 
the Chi-square error (CSE), defined in Eq. (26) is used as the 
second measure just as Liu et al. did in [34]: 


n oy. AR 
ceay ee (26) 


{Ti Xa 


As featured in Eq. (26), CSE is actually a relative error between 
(Xoi—Xci)? and xali=1,2,...,n), thus here the third criterion which 
can be used to measure the relative error between (xoi— Xci)? and 
Xoi (i=1,2,...,n) is adopted. Because of its similarity with the CSE, 
it has been called the similar Chi-square error (SCSE): 


n 2 
(Xoi—Xci) 
SCSE = a AR 27 
In addition, the root mean square error (RMSE) is also 
employed to evaluate the performance of different distributions: 


1/2 
— 1 í P . 2 
RMSE = f XO oiai) | : (28) 


i=1 


7.3. Simulation and comparison results 


To analyze the wind speed distribution in the area described in 
Section 3 by the three distributions, i.e., Weibull, Logistic and 
Lognormal, the following procedures are employed in dividing 
wind speed ranges: 


Step 1: Marking the minimum and the maximum wind speed 
data as MI and MA. 

Step 2: Generating several intervals with the equal length 
1 like this: [Floor(MI),Floor(MI)+ 1),[Floor(MI)+ 1,Floor(MI)+ 2), 


... [Ceil(MA)—1,Ceil(MA)), where Ceil(MA) is the nearest inte- 
ger larger than MA and Floor(MI) is defined as itself if MI is an 
integer; otherwise, it is set as the nearest integer smaller 
than MI. 


By the above range division method, the yearly wind speed 
data are divided into different intervals and the number of wind 
speed data located in each interval is calculated. Fig. 4 gives the 
frequency of the wind speed located in each wind speed interval 
for the twelve months throughout the year. 

From Fig. 4(a) it can be observed that 16 intervals are 
generated in all: [2,3),[3,4), ...[17,18). In addition, the wind speed 
data located in the interval [8,9) are the greatest, and in the 
interval [9,10) are the second greatest, while the wind speed data 
fall in range [17,18) are the least. Fig. 4(b) is the x-y view of 
Fig. 4(a). It shows that all of the 16 intervals contain the wind 
speed data in December. A y-z view of Fig. 4(a) is shown in 
Fig. 4(d). It can be observed that the maximum frequency of wind 
speed located in interval [8,9) is about 260, while the minimum in 
the range [17,18) is close to 0. 

By the statistical number given in Fig. 4, the probability density 
function and the cumulative distribution function in each month can 
be easily obtained. Now applying the three different distributions 
mentioned in Section 6.1 to wind speed modeling, the pdf and cdf of 
these three distributions are obtained. Fig. 5 shows the monthly pdf 
and cdf of the observed wind speed and three distributions in four 
different months: February (in Winter), May (in Spring), August (in 
Summer) and November (in Autumn). As presented in Fig. 5(a)-(d), 
the probability for wind speed data located in ranges [10,11), [8,9), 
[8,9) and [9,10) are the largest for February, May, August and 
November, respectively. The maximum probability values calculated 
by the Logistic pdf for these months are all much closer to the 
observed probability values than those calculated by the other two 
distributions. For pdf lines in May displayed in Fig. 5(b), the observed 
and the calculated maximum probability values by the Logistic pdf 
are greater than 0.2, while for the other three months, the maximum 
probability values are smaller than 0.2. Additionally, the maximum 
wind speed in May is located in the interval [11,12), while the 
maximum wind speed data for the other three months are all larger, 
and the one in February which is located in the range [17,18) is the 
largest. In addition, the cdf lines plotted by the three distributions for 
February fit the actual cdf line better than for other months. 

Apart from the monthly wind speed distribution modeling, the 
yearly wind speed distribution analysis is also quite meaningful. 
Thus Fig. 6 provides the pdf and cdf comparison between the 
observed wind speed data and data calculated from the three 
distributions. Similar to May and August, the probability for wind 
speed data located in the interval [8,9) is the maximum. 

Since figures only provide a rough view to different models or 
distributions, thus the numerical calculations for the three dis- 
tributions are provided to give a clear description. Error values 
between the observed and the calculated values for the 12 
months and the whole year are listed in Table 5. As seen, there 
are several incalculable values when the MRE or SCSE is selected 
to be the performance evaluation criterion. This is because once 
zero value appears in the sequence {x,;}/_ ,, the definition of the 
MRE or SCSE makes no sense. Thus these two error criteria have 
drawbacks. However, for CSE, the same defect also exists when 
the zero value appears in the sequence {x,;}7_ ,, while the fourth 
evaluation criterion (RMSE) can be put to use in any cases. 

Results in Table 5 reveal that the MRE obtained by the Weibull 
distribution is the smallest for May and July as compared to the 
other two distributions. Apart from February, April and Decem- 
ber, in which the smallest MRE values are obtained by the Logistic 
distribution, for other months and the whole year, the Lognormal 
distribution performs best under this error evaluation criterion. 
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When CSE, SCSE or RMSE is selected to be the error evaluation 
criterion, the situation is similar: For some months, the Weibull 
distribution performs best in fitting the wind speed distribution, 
while for other months, the best distribution is the Logistic or the 
Lognormal distribution. 

As a result, under different evaluation criteria and months, the 
three distributions perform discriminately. In Table 5, the three 
distributions are compared 13 times (12 months and the whole 
year) under each judgment criterion, thus for each time, three error 
values are obtained. Recording the minimum and the maximum 
value as ‘S’ and ‘Ll’, and the remaining one as ‘M’, then the 
performance frequency of the three distributions under each eva- 
luation criterion is calculated, the results are plotted in Fig. 7. As 
shown in the MRE of Fig. 7, the best performance frequency for 
Weibull, Logistic and Lognormal distributions are 2, 6 and 3 respec- 
tively. These values change to 5, 5 and 3 for CSE, 3, 5, 3 for SCSE and 
3, 4, 6 for RMSE. Consequently, as compared to the Logistic and 
Lognormal distributions, the performance of the Weibull distribu- 
tion is not as good as the one mentioned in [34], which used the 
Weibull function to analyze the wind speed distribution finally for 
its best ability. Thus, there is no consistent conclusion to tell which 
function is the best in modeling wind speed distribution. Function 
performance should be analyzed and understood under specific 
conditions. Optimal function should be determined according to the 
time, area resource of the sampled data or the evaluation criterion 
which has been taken into account preferentially. For example, if 
one prioritizes the RMSE, the Lognormal function will be the 
optimal; while if only SCSE and RMSE are considered and they 
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occupy the same weight, the frequency of the appearance of ‘S’ 
appears for the three distributions should be computed firstly. The 
calculated values are 6 (3+3=6), 9 (5+4=9) and 9 (3+6=9) 
respectively, thus both of the Logistic and Lognormal distributions 
can be regarded as the best. 

In this paper, considering the drawback of the MRE and the 
SCSE criterions in calculating the errors, these two criterions have 
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Fig. 6. Yearly pdf and cdf of the observed wind speed and three distributions. 
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Fig. 5. Monthly pdf and cdf of the observed wind speed and three distributions in four different months: (a) February, (b) May, (c) August, (d) November. 


J. Wu et al. / Renewable and Sustainable Energy Reviews 21 (2013) 215-228 


been given lower weight coefficient than the other two criterions 
CSE and RMSE. The weight coefficients for MRE, CSE, SCSE and 
RMSE are 1/6, 1/3, 1/6 and 1/3 respectively. So the final perfor- 
mance of the three distributions is 


Weibull: 1/6*2+1/3*5+4+1/6*341/3*3=7/2; 

Logistic: 1/6*6+1/3+*5+4+1/6+5+41/3+4=29/6; 

Lognormal: 1/6*3+1/3*341/6*3+1/3+*6=4. 

For 29/6 > 7/2 and 29/6 > 4, thus the Logistic distribution is 
the best distribution in fitting the yearly wind speed in this area, 
and the following analyses are all based on this best distribution. 


8. Wind energy calculation and analysis 
8.1. Wind energy calculation 
The most probable wind speed for a particular location with the 


given probability function can be computed according to the calcula- 
tion of the maximum value of the Logistic pdf by of(v)/ov=0, by 


Table 5 
Comparison results of three distributions. 
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which the following formula can be obtained: 
epiw -Dderw L o, 29) 
0° -{1+exp[—(v—1) /d} 
thus the most probable wind speed can be expressed as 
Uup = V. (30) 


The cumulative Logistic distribution function enables the 
calculation of the probability of the wind speed exceeding the 
value u: 


F(v > u) = 1—F(v < u) 
1 : 1 
1+exp[—(v—v)/d] 1+exp(T/ð)` 


(31) 


For a wind turbine designed with a given cut-in speed v; and 
cut-off speed vo, the operating probability of the wind turbine, 
named here availability factor, is equal to the probability of wind 
speed between v; and vo, calculated by 


Ar = F(V; < V < Vo) = F(v < Vo) —F(V < vi) 


Period MRE CSE SCSE RMSE 
Weibull Logistic Lognormal Weibull Logistic Lognormal Weibull Logistic Lognormal Weibull Logistic Lognormal 
January 0.3995 0.3353 0.6762 0.0627 0.0652 0.3609 0.0830 0.0565 0.2446 0.0190 0.0170 0.0395 
February 0.2803 0.2881 0.2459 0.0952 0.0973 0.0667 0.0949 0.0977 0.0592 0.0242 0.0231 0.0200 
March 0.3869 0.2197 0.2518 0.0670 0.0319 0.2075 0.0706 0.0309 0.0246 0.0229 0.0156 0.0113 
April 1.0433 0.7923 0.7118 0.4016 0.2480 0.1922 0.5142 0.3037 0.2928 0.0384 0.0369 0.0363 
May 0.2134 0.2512 0.3661 0.0535 0.0871 0.1774 0.0606 0.0888 0.1726 0.0272 0.0334 0.0434 
June =b =i =» 3.1385 0.1324 0.3173 =» -» => 0.0302 0.0237 0.0144 
July 0.1369 0.2040 0.2432 0.0252 0.0538 0.0692 0.0226 0.0487 0.0539 0.0178 0.0251 0.0228 
August 0.3043 0.2850 0.3932 0.0767 0.0959 0.2118 0.0837 0.0702 0.0980 0.0218 0.0216 0.0211 
September 0.4460 0.3611 0.3768 0.1738 0.1129 0.1719 0.1414 0.0858 0.1395 0.0270 0.0229 0.0313 
October = -- =b 0.1246 0.1326 0.0659 =b -P an 0.0271 0.0307 0.0180 
November 0.2210 0.2079 0.2772 0.0372 0.0821 0.0969 0.0375 0.0639 0.0822 0.0149 0.0205 0.0244 
December 0.5190 0.3746 0.3376 0.1070 0.0699 0.2020 0.1276 0.0745 0.0861 0.0234 0.0198 0.0275 
Yearly 0.5227 0.4192 0.6245 0.0693 0.0294 0.0520 0.0552 0.0373 0.0454 0.0130 0.0103 0.0112 
P Incalculable value. 
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Fig. 7. Frequency of the three distributions under each judgment criterion. 
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1 1 
~ 1+exp[(—v,—v)/d] 1+ exp[—(1;—1)/6] 


(32) 


8.2. Capacity factor 


The capacity factor is a very significant index of productivity of 
a wind turbine. It is the ratio of the actual energy produced in a 
given period, to the hypothetical maximum possible, i.e. running 
full time at rated power. In this paper, the capacity factor is 
calculated by the following expression: 


Ew 

ER" (33) 

where 

Ew=T | Pre vifv) do, a 

0 

Pr, Vr <V < Vo, 
v—vi 

Pwrg(V) =? P, T; , Vi <V < Vr, (35) 
0, v<vjorv>vVp, 

Er =TP;. a6 


Vi, Vp, Vo, P, and T represents the cut-in, rated, cut-off speed, the 
rated power and the time. In this paper, f(v) is the probability 
density function of the Logistic function. 


8.3. The turbine efficiency 


The turbine efficiency 7, which is defined as the ratio of the 
recoverable energy on the aerogenerator and the available energy 
at Betz limit, can be simply calculated by the following formula 
[37]: 


(37) 


where t the time period, A is the blade sweep area of wind 
turbine, P, is the rated power of the wind turbine and P; is defined 
as 


t 1_3/h\* 
Ph = T0002? ¥m (xs) , (38) 
in which Vm is the mean wind speed of this time period at the 
standard height 10m, h is the hub height, « is the roughness 
factor similar as Eq. (1), and the corrected air density p (kg/m?) is 
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8.4. Wind energy results and analysis 


Wind and wind turbine characteristics for different time 
periods have been summarized in Table 6. As seen, all of the 
availability factors in this area are larger than 0.98. The capacity 
factor presents a similar trend with the mean wind speed for each 
time period: the capacity factor in February is the largest. Mean- 
while, the mean wind speed in this month (see Table 1) is also the 
maximum. Similarly, the smallest capacity factor value and the 
minimum most probable speed appear in the same month (May). 
So in the month with the larger mean wind speed, the wind 
turbine could utilize more wind energy. Furthermore, the smal- 
lest value for the turbine efficiency is 0.6555. According to Betz, 
under the assumption that the turbine swirl and transmission 
losses were ignored, the theoretically maximum power that can 
be extracted from the wind is 59% of the wind power available in 
the wind. Therefore, for any wind turbine, the wind turbine 
efficiency should not exceed 0.59 [37]. Thus, the wind energy 
potential in this site of Inner Mongolia appears to be abundant. 

To analyze the wind energy output, six other different types of 
wind turbine generators are adopted. The specifications of these 
turbines are given in Table 7, among which the maximum cut-in, 
rated and cut-off speeds are 5, 15 and 27 m/s, respectively. The 
maximum rated power is as high as 3000 kW. Meanwhile, two 
wind turbines (CONE-450 and YT/850) with the same hub height 
(50 m) with the given wind turbine are selected to survey the 
influence of different rated speed to the availability and capacity 
factors. 

The availability and capacity factors for different wind tur- 
bines are listed in Table 8. As seen, the availability factors of the 
S43/600, CONE-450 and YT/850 are the same for all of the months 
and the whole year. This is because the cut-in and cut-off of these 
three wind turbines are the same as well as the hub height. It can 
be concluded that the availability factor is only relative to the cut- 
in and cut-off speed, but has no relation with the rated speed and 
the rated power. In addition, for wind turbines with the same cut- 
in and cut-off speed, the higher the hub height, the larger the 
availability factor. This can be viewed by the availability factor 
values of the NODEX-N 70/1500 and SL 3000/90 wind turbines. 
Capacity factor values for different wind turbines in Table 8 


Table 7 
Specifications of different wind turbines. 


determined by Number Type P/kW H/m v;/ (m/s) v,/ (m/s) vo/ (m/s) 
P 1 CONE-450 450 50 3 12 25 
P= = (39) 2 NM600/43 600 56 4 15 25 
Ra*T 3 YT/850 850 50 3 15 25 
= Pirie Sek a 4 NODEX-N70/1500 1500 70 3 13 25 
where P and T are the mean air pressure and the mean air 5 GAMMA-60 1500 66 5 133 27 
temperature in the time period t with the unit of N/m~ and Kelvin, 6 SL 3000/90 3000 90 3 13 25 
respectively and Rg=287. 
Table 6 
Wind and wind turbine characteristics for different time periods. 
Item Period 
January February March April May June July August September October November December Yearly 
Logistic scale parameter (m/s) 1.0863 1.4730 1.0184 1.1528 1.0399 1.0504 1.1165 1.2996 1.3165 1.3791 1.4659 1.1522 1.3087 
Most probable speed (m/s) 8.9439 10.2558 8.6957 8.2298 7.4566 7.5508 7.8993 9.6786 9.5821 9.5605 9.6972 8.7084 8.8444 
Availability factor 0.9958 0.9927 0.9963 0.9894 0.9864 0.9870 0.9877 0.9942 0.9933 0.9915 0.9897 0.9930 0.9886 
Capacity factor 0.4954 0.6006 0.4748 0.4366 0.3725 0.3803 0.4093 0.5554 0.5475 0.5455 0.5561 0.4760 0.4873 
The corrected air density 1.1647 1.1686 1.1692 1.1772 1.1834 1.1898 1.1999 1.1975 1.1901 1.1824 1.1760 1.1695 1.1807 
The turbine efficiency 0.9917 0.6555 1.0748 1.2594 1.6843 1.6133 1.3971 0.7611 0.7892 0.7998 0.7706 1.0699 1.0116 
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Table 8 
Availability and capacity factors for different wind turbines. 
Factor Number Period 
January February March April May June July August September October November December Yearly 
Availability factor 1 0.9958 0.9928 0.9963 0.9894 0.9864 0.9870 0.9877 0.9942 0.9933 0.9915 0.9897 0.9930 0.9886 
2 0.9901 0.9864 0.9907 0.9764 0.9672 0.9689 0.9720 0.9881 0.9864 0.9833 0.9807 0.9843 0.9770 
3 0.9958 0.9928 0.9963 0.9894 0.9864 0.9870 0.9877 0.9942 0.9933 0.9915 0.9897 0.9930 0.9886 
4 0.9963 0.9933 0.9968 0.9906 0.9881 0.9886 0.9891 0.9947 0.9939 0.9923 0.9906 0.9938 0.9897 
5 0.9782 0.9758 0.9785 0.9511 0.9273 0.9315 0.9409 0.9769 0.9740 0.9691 0.9656 0.9672 0.9562 
6 0.9966 0.9937 0.9971 0.9914 0.9892 0.9896 0.9901 0.9951 0.9944 0.9928 0.9911 0.9943 0.9905 
Capacity factor 1 0.6539 0.7636 0.6290 0.5777 0.4953 0.5055 0.5428 0.7205 0.7107 0.7061 0.7150 0.6280 0.6385 
2 0.4630 0.5792 0.4402 0.3987 0.3282 0.3367 0.3687 0.5294 0.5207 0.5186 0.5304 0.4418 0.4548 
3 0.4954 0.6007 0.4748 0.4366 0.3725 0.3803 0.4093 0.5554 0.5475 0.5455 0.5561 0.4760 0.4873 
4 0.6333 0.7436 0.6091 0.5606 0.4820 0.4917 0.5273 0.6994 0.6899 0.6859 0.6953 0.6088 0.6197 
5 0.5170 0.6564 0.4873 0.4329 0.3407 0.3518 0.3938 0.5995 0.5884 0.5847 0.5977 0.4891 0.5049 
6 0.6645 0.7710 0.6401 0.5899 0.5094 0.5193 0.5558 0.7292 0.7196 0.7151 0.7236 0.6391 0.6492 
illustrate that for each month and the whole year, the capacity Acknowledgment 


factor of the SL 3000/90 is the largest, and the one of the CONE- 
450 is the second largest. The reason is that the capacity factor is 
influenced by the cut-in, rated and cut-off speeds, as well as the 
probability density function. This can be clearly seen by the 
simplification type of Eq. (33) shown as below: 


* y—v, 
Crs I t F(u)du+ Ar. (40) 
vi i 


Vr— 


9. Conclusions 


In this paper, wind resource in one of the stations in Inner 
Mongolia, China, is analyzed. First, the wind speed probability 
distribution is modeled by three probability density functions: the 
Weibull, Logistic and Lognormal functions, among which the shape 
parameter in the Weibull function is estimated through three PSO 
algorithms and 18 DE approaches, and the one with the smallest 
minimum iteration number under the same fitness level is selected to 
determine the final Weibull shape parameter. Then four error 
evaluation criteria are adopted to choose the best probability density 
function by giving them different weight, and the best Logistic 
function is used in the consequent wind energy assessment. Capacity 
and availability factors are two necessary aspects in wind energy 
assessment. Thus, the two factors of six other different wind turbines 
are compared with the selected wind turbine. The conclusions are 
summarized as below: 


(a) There is no function that performs best under any conditions 
in wind speed distribution modeling. So function performance 
should be analyzed and understood under specific conditions. 
Optimal function should be determined according to the time, 
area resource of the sampled data or the evaluation criterion 
judgment criterion which is taken into account preferentially. 

(b) Performance of the intelligent parameter optimization meth- 
ods, such as the PSO and DE, is not always better than that of 
the conventional moment estimation approach in searching 
the optimal parameter values. 

(c) The available wind speeds in this area for each month and the 
whole year are larger than 98%. 

(d) The capacity factor presents a similar trend with the mean 
wind speed for each time period, thus in the month with the 
larger mean wind speed, the wind turbine can utilize more 
wind energy. 

(e) This area is suitable for a wind farm foundation. 
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this paper. 
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